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Abstract. In the Randall-Sundrum model with one brane, we found the 
approximate and exact solutions for gravitational potentials and accelerations of 
test bodies in these potentials for different geometrical configurations. We applied 
these formulas for calculation of the gravitational interaction between two spheres 
and found the approximate and exact expressions for the relative force corrections 
to the Newton's gravitational force. We demonstrated that the difference between 
relative force corrections for the approximate and exact cases increases with the 
parameter I (for the fixed distance r between centers of the spheres) . On the other 
hand, this difference increases with decreasing of the distance between the centers 
of the spheres (for the fixed curvature scale parameter I). We got the upper limit 
for the curvature scale parameter I < 10 ^m. For these values of I, the difference 
between the approximate and exact solutions is negligible. 



PACS numbers: 04.50.-h, 11.25.Mj, 98.80.-k 
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1. Introduction 

The idea of the multidimensionahty of our Universe has been attracting continuous 
interest for many years. It takes its origin from the pioneering papers by Th. Kaluza 
and O. Klein [T], and now the most self-consistent modern theories of unification such 
as superstrings, supergravity and M-theory are constructed in spacetime with extra 
dimensions [2]. In Kaluza-Klein models, our spacetime is effectively four-dimensional 
due to compactness and smallness of the extra dimensions (internal spaces). The size 
of the extra dimensions is restricted by the electroweak scales 10~^^ cm. However, 
our spacetime can be effectively four-dimensional even in the case of infinite extra 
dimensions. This interesting scenario is realized in recently proposed brane world 
models (see, e.g., the reviews [3j|4]). Here, matter fields from the Standard Model are 
trapped to a three-dimensional submanifold (brane) embedded in the fundamental 
multidimensional space (bulk), but gravity may move in the bulk. Localization of 
massless gravitons on a brane results in effective four-dimensional Einstein gravity in 
the low energy limit. Certainly, large and infinite extra dimensions are potentially 
detectable. This was one of the main reasons for the great interest in this scenario. 
Therefore, it is very important to suggest experiments which can reveal such extra 
dimensions. 

In our paper, we consider the scenario that was first proposed in 5 . Here, 
the brane is embedded in the five-dimensional anti-DeSitter spacetime, which allows 
the extra dimension to be infinite. A negative bulk cosmological constant A5 and 
a brane tension a are fine tuned to each other. Clearly, this is a very simplified 
scenario. However, it gives a possibility to reveal some general features of the brane 
world models, in particular, the localization of the massless graviton on the brane 
that restores the Newtonian limit on the brane at large distances from the gravitating 
matter source. It was shown [5^ that at distances greater than a curvature scale 
of anti-DcSittcr spacetime r ^ I ^ jAsj"^/^, the gravitational potential takes an 
approximate form with a cubic additive ^ to the usual Newtonian potential 

1/r. This approximate solution is much simpler than the exact one that makes 
the investigation of the effects of the extra dimension much easier. In some papers 
(see, e.g., [Bl[7|) this approximation was used to calculate the gravitational interaction 
between gravitating test bodies of different geometrical form. But we should analyze 
the difference between the approximate and exact solutions to find out where the 
application of the approximate solution is appropriate. This is one of the main 
motivations of this work. To perform such analysis, we obtain two types of solutions 
(approximate and exact) for gravitational potentials and accelerations of test bodies 
in these potentials for different geometrical configurations. Then, we apply these 
formulas to the most interesting for experiments case of gravitational interaction 
between two massive spheres. We calculate approximate and exact corrections to 
the Newton's gravitational force and show that the difference between relative force 
corrections for the approximate and exact cases increases with the parameter I (for 
the fixed distance r between centers of the spheres) . On the other hand, this difference 
increases with decreasing of the distance between the centers of the spheres (for the 
fixed curvature scale parameter /). The relative force corrections also allow us to get 
the experimental constraint on the curvature scale parameter: / < 10 /im. To get it, we 
use the results of the table-top inverse square law experiments for the measurements 
of the Newton's gravitational constant. This is one of the main results of our paper. 
The paper is structured as follows. In section 2 we describe briefiy the Randall- 
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Sundrum model with one brane. Here, we consider non-relativistic limit of this model 
and present approximate and exact solutions for the gravitational potential on the 
brane. These formulas are applied to some practical problems in section 3 to get 
approximate and exact expressions for the gravitational potential and acceleration of a 
point mass for these problems. In section 4 we investigate the gravitational interaction 
of two spherical shells. Then, in section 5 we compare the relative corrections to the 
gravitational force between two spheres in approximate and exact cases. Here, we also 
get the constraint on the curvature scale parameter in the Randall-Sundrum model. 
A brief discussions of the obtained results is presented in the concluding section 6. 



2. Non-relativistic limit of Randall-Sundrum model 

The one-brane Randall-Sundrum metrics is [5] 



ds^ = exp y p j r^^^^dxt'dx" - d^^ , (2.1) 

where 77^1/ is the flat four-dimensional spacetime metrics and the parameter I is defined 
via the 5-dimensional cosmological constant: 

A5 = -^, (2.2) 

i.e. I is the curvature scale of 5-dimensional anti-DeSitter spacetime. The brane is 
embedded in this spacetime at ^ = and has fine tuned tension 

where G5 is the 5-dimensional gravitational constant. In one-brane Randall-Sundrum 
scenario, the extra dimensions is infinite: ^ G (— oo,-|-c»). 

Now, we want to probe this model with the help of the gravitational terrestrial 
experiments, e.g., the inverse square law experiments. Certainly, this is the case of 
non-relativistic limit of the model. In this limit, we need to get the gravitational 
potential ip{r) on the brane. Following, e.g., the calculations in [4], we obtain 

00 

ip[r) ; / dm^p^^[l) exp[~mr) , (2.4) 

rl r J 



where r is the magnitude of the radius vector on the brane and 



rn?y/^ Yi{ml)J2{ml) - Ji{ml)Y2{rhl) 



where J and Y are Bessel functions of the first and second kinds, respectively. In the 
short and long distance limits the equation (|2.4p reads, respectively: 



( \ ^ G5TO 1 Gsm 

Vin w « , r < Z (2.6) 

ri r irr irr 



and 

where we have introduced the Newton's gravitational constant 



V{r)- —[^ + -) ' (2-7) 

r \ r J 



Gn — -y- (2.' 
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and the parameteid 

P 

a - - . (2.9) 

Obviously, \2.Q\ corresponds to the strong deviation from the Newtonian gravity 
but the formula (|2.7p describes the smooth transition to the Newtonian limit. The 
exact expression p.4|) (the solid line) and its asymptotes ()2.6|) and (|2.7|) (the short- 
dashed and long-dashed lines, respectively) are depicted on figure 1. Here, we 
introduce the dimensionless distance argument rj ~ r/l and dimensionless potentials 
¥>iv) = v{r)/{GNm/l). 

3. Applications 

Now, we want to apply the obtained formulas to terrestrial gravitational experiments. 
Obviously, the gravitational field on the Earth should not considerably differ from 
the Newtonian one. Therefore, we should use either the exact expression (|2.4[) or the 
approximate formula (|2.7p . Therefore, we shall get two classes of solutions: exact and 
approximate, respectively. Obviously, the approximate formula (|2.7p looks much more 
simple. However, it is necessary to check the deviation of expressions based on it from 
the exact ones for real gravitational experiments. This is one of the main aims of the 
paper. 

3.1. Infinitesimally thin shell 

Let us consider first an infinitesimally thin shell of the mass m = AnR^a, where R 
and a are the radius and the surface mass density of the shell. Then, the gravitational 
potential of this shell in a point with the radius vector r (from the center of the shell) 

f It is worth noting that in the pioneering paper [5] a = Z'^. The brane-bending effect 8 gives 
2/2/3. In the paper 9], the authors have pointed out that different schemes of regularization 
result in different values of a. In our paper we follow calculations in 4, where a = /2. 
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for the approximate solution is 



and 



Lp{r > R) 



(p{r < R) 



R 



i?2 



(3.1) 

(3.2) 



Obviously, these expressions are divergent when r ^ R: (p(r ^ R) 
—GNunal [2i?^|r — — > —oo. In the case of the exact solution we get 



ip{r > R) 
and 

ip{r < R) 



2R 



dm 



Gnttt- 
R 



1 + — / dm 

2r 



Yi{ml)J2{ml) - Ji{ml)Y2{ml) 
{Jl{ml) + Yl{ml)Yl^ 



Yi{ml)J2{ml) - .h{ml)Y2{ml) 
{Jl{ml) + Y^{ml)Y/-^ 



sinh(mi?)e 



(3.3) 



sinh(mr)e 



— rriR 



(3.4) 

It is not difficult to verify that these exact solutions are also divergent when r R. 

Formulas p.2p and (|3.4p demonstrate that inside of the shell the gravitational 
potential is not a constant. Thus, a test body undergoes an acceleration (see (|3.6p 
and p.Sp below) in contract to the Newtonian case, i.e. the Birkhoff's theorem is 
violated. The acceleration outside and inside of the shell is 



'-fir>R) 
dr 



G^m 



and 



1 + a 



2ar 



3r' - 



(r2 - i?2)2 



d^, GNm 

dr^'^ ' R {R^-r^Y 



which is divergent when r ^ R: -^{r ^ R) ^ "^iSfp^ 
lower signs correspond to r > R and r < R, respectively). 
In the case of exact solutions we have 
Gatto 



(3.5) 
(3.6) 

=pcx) (the upper and 
(3.7) 



GnitiiI 



2Rr^ 



dm 



Yi{ml)J2{ml) - Ji{ml)Y2{ml) 
{Jl{ml) + Y^{ml)Y/^ 



sinh(mi?)(l + mr)e 



and 



'f[r<R) 
dr 



GNml 
2Rr^ 



dm 



Yi(ml)J2{ml) - Ji{ihl)Y2{fhl) 
{Jl{ml) + Y^{ml)Y/-^ 



(mr cosh(TOr) — sinh(mr))£ 



(3.8) 



-rhR 



The exact solutions (j3.7p and (j3.8p are also divergent for r 
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3.2. Spherical shell of finite thickness 

Here, we consider a spherical shell of the inner radius Ri and the outer radius R2 
and the mass m — (47rp/3) (i?2 — with a constant volume density ^ For this 
geometry, the approximate gravitational potential reads 



(p(r > i?2) = - 
and 



Gntti 



3a 



i?.-i?, + ^inj!:±^#-^ 

- 2 {r-R2){r + Ri) 



ifiir < Ri) = -2ttGnP [ R2 - Ri 



Rt_ 

Rl- 



(3.9) 



(3.10) 



These expressions are logarithmically divergent in the vicinity of Ri and R2: ip{r > 



R2) 



ln(r — i?2) ^ —00 for r — > i?2 and (p{r < i?i) — > 2tiGnPol ln(i?i — r) 



2{Rl-Rl) 

00 for r — > _Ri . 

For the exact solution we get 



(p(r > i?2) = - 
3 

and 



1 



dm 



Yi(ml)J2{ml) - .h(ml)Y2{ihl) 



n 2 



{Jl{ml) + Y^{ml)YI-^ 

3 p3\ cosh(TOi?) — sinh(mi?)] 
(it2 — Ri) \ 



<p(r < i?i) = -AttGnP 




Yi{ml)J2{ml) - Ji{ml)Y2{ml) 
{Jl{rhl) + Y^{ml)Yl^ 



(3.11) 



■ sinh(mr) 



X {fnR + l)e '' 



I-R2 



(3.12) 



In contrast to the approximate formulas p.9p and p.lOp . these exact expressions are 
convergent in the limits r — J> Ri, R2. 

The acceleration of a test body outside and inside of the shell is 



dip 



dr 
and 

dip 
dr 



{r > R2) 



Gi^m 



1 



3q! 



i?| - Ri 



2{R2-Ri){r'+RiR2) 

Ki — K2 + r 



{r^ - Rl){r^ - R. 



(3.13) 



{r < Ri) = AnGNpar—j 
[R2 



Rl 



Rl 



r^){Ri~r^) 



These formulas are divergent in the limits r 



-> -00 for r i?2 and ~^{r < Ri 



2{Rl-Rl){r-R,_) 

For exact solutions we obtain 



(3.14) 

Ri,R2- -^{r > R2) ^ 
^ +00 for r^Ri. 



^lr>U, 



G 



GnuiP 
2r2 



dm 



Yi{ml)J2{ml) - Ji{ml)Y2{ml) 



X -ire-™'- (1 + mr) 



{J^{ml) + Y^{rhl)y/^ 
R2 



—3 —3- [_Rto cosh(mi?) — sinh(?7ii?)] 

(R2 ~ Rl) ^ 



(3.15) 



§ It is clear that the hmit R2 — > Ri is incorrect because, for constant p, it results in vanishing m 
and, vice versa, for fixed m the volume density p goes to infinity. Therefore, such a naive limit does 
not provide us the correct transition to the formulas from the previous subsection. 
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and 



dtp 
dr 



{r < Ri) 



dm 



Yi{ml)J2iml) - Ji{ml)Y2{ml) 



{Jl{ml) + Y^{ml)Yl^ 



X -r-j- (mr cosh(mr) — sinh(mr))(TOi? + l)e ™ | 



-rnRiB.2 



These integrals are divergent in the vicinity of i?i and i?2 ■ 



(3.16) 



3.3. Sphere 

Obviously, all formulas for a sphere of the radius R and the mass m = 47rpi?^/3 with 
a constant volume density p can be easily obtained from the equations p.9p . p.lip . 
p.l3p and p.lSp with the help of the evident substitutions: Ri = and i?2 = ^• 



4. Gravitational interaction of two spherical shells 



Let us consider now two spherical shells with radii R2 > Ri and the mass m ~ 
(477/3/3) (i?2 ~ Rf) '^'"^t shell and radii R2 > R'l and the mass m' — 

(47ryo'/3) (^R2 — Ri^ for the second shell. Then, the potential energy of gravitational 
interaction between these shells for the approximate solution reads 

, GNfnm' 2'K^GNpp'a 

U{r) = - 



r r 
„2 fW I o\2 



1 4 I 1 2 / d'2 



R' 



, r- ~ lit' + ur 2 
X In — — -frr + -r 



- -r^R'R 
3 

where r ^ i?2 



(i?^ + RY 
{R'-RY ' 3 

R'^R - R'R^ 



■^.,Jr + RY-R"- 



(r - i?)2 - i?'2 

\r=Ri\r'=r' ' 



i?'' In 



i?2 



(r - i?')2 - 



(4.1) 



/(i?,i?')lfl=Sli?'=ffi =/(-fi^2,i?2 



is the distance between the centers of the shells and 



U{r) 



In the case of the exact solution we get 
GNfnm' GNmm' 9 



dm 



[Rl - RlYR'2' ~ R[^) 



Yi{ml)J2{ml) - Ji{ml)Y2{ml) 



{J^{ml) + Y^{ml)y/^ 



X [mRcosh{mR) — sinh(m_R)] • [mi?cosh(mi?) - sinh(mi?)] . (4.2) 

The additional analysis shows that both of these expressions (|4.ip and (|4.2I) are 
convergent in the limit r — >■ i?2 + -R2 • 

With the help of these formulas, we can obtain the absolute value of the 
gravitational force between two shells: 
dU GNmm'^ 
dr 7 



F(r) 



(4.3) 



where 5f defines the relative deviation from the Newtonian expression Gf^mmJjr'^. 
For the approximate and exact solutions we have, respectively: 



5f= - 



9a 



8(i?3-i?3) (i?^3_^/3) 



In- 



{R! + RY 
{R! - RY 



1 4 

r 

4 



R' 
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and 

5f^ 



9P 



r^R'R + R''R + R'R''\\lzl\\irJf 



(4.4) 



2 {Ri ~ Ri) {R:i ~ 



dm 



Yi(ml)J2{ml) - Ji(ml)Y2{ml) 



{Jf{ml) + Y^^{ml)y/-^ 



X (1 + mr) [mRcosh{mR) - sinh(mi?)] • [rhR' cosh{mR') - sinh(mi?')] 1^? • 

'(4.5) 

These relative corrections 5f are also convergent in the limit r — i?2 + ^2- 
the limit of large separation between the shells r ^ Ri,2,R[ 2 '^^ obtain from (|4.4p 
5f = 3a/r2 +0(l/r3). 



5. Constraints 

The obtained above formulas can be used for the experimental restrictions on the 
parameters of the model. In our case, it is the curvature scale I. To get it, we can use 
the inverse square law experiments for two spheres. The potential energy of interaction 
and the gravitational force between two spheres follow from the previous section with 
the help of the substitutions: Ri = R[ = and R2 = R,R'2 = R' ■ For example, the 
relative corrections to the gravitational force in approximate and exact cases read, 
respectively: 



Sf= - 



9a 



In- 



{R' + Rf 



8R^R'^ 
^R'R + R'^R 



■-{R' - 
R'R'} 



i?)2 



4 



2 ^ 



(5.1) 



and 



9f 



2i?3i?'3 



dm 



Yi{ml)J2iml) - .h(ml)Y2{ml) 



{Jl{ml) + Y^{ml)Y/^ 
sinh(mi?)] • [mi?' cosh(mi?') — sinh(mi?') 



(1 + mr) 



(5.2) 



X [mi?cosh(mi?) 

We remind that a — /^/2 (see the equation (j2.9|) '). For definiteness, we shall use 
the parameters of the spheres from the Moscow Cavendish- type experiment |10) : 
Ri ~ 0.087 cm for a platinum ball with the mass mi = 59.25 x 10"^ ~ 0.206 cm 

for a tungsten ball with the mass m2 = 706 x 10"^ g and the distance between their 
centers r = 0.3773 cm. 

It is clear that the use of the approximate solution for the gravitational interaction 
force makes the calculations much easier. But we should analyze the distinction 
between the approximate and exact solutions to find out where the application of 
the approximate solution is appropriate. The difference between the approximate and 
exact solutions for the relative force corrections is shown on figures 2 and 3. These 
figures demonstrate that the difference between these solutions increases with the 
parameter / (for the fixed distance r between centers of the spheres) (figure 2). On 
the other hand, this difference increases with decreasing of the distance between the 
centers of the spheres (for the fixed curvature scale parameter /) (figure 3). 

Now, we want to estimate the curvature scale parameter I with the help of 
our formulas (|5.ip and (|5.2|) . To get it, we can use the value of the Newton's 
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Figure 2. Relative gravitational force corrections 115.21 1 (the solid line) and 115.111 
(the dashed line) as functions of the curvature scale parameter I in the case of the 
distance between the centers of the spheres r = 0.3773 cm. 
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Figure 3. Relative gravitational force corrections 115.21 1 (the solid line) and 115.111 
(the dashed line) as functions of the distance between the centers of the spheres 
in the case I = 10~^ cm. 

gravitational constant Gn- As it follows from figure 2 in the 'CODATA Recommended 
Values of the Fundamental Constants: 2006', the most precise values of Gn were 
obtained in the University Washington and the University Ziirich experiments [TTl [T^ . 
They are GAr/lO-^m^kg^^s-^ ^ 6.674 215 ± 0.000 092 and 6.674 252 ± 0.000 124, 
respectively. The relative errors AGn/Gn show the accuracy of the measurements 
of the gravitational constant in the inverse square law experiments. If the correction 
6f due to the extra dimension is greater than these values, then we can detect the 
deviation from the Newton's law. Up to now, there is no experimental evidence for such 
deviations. Therefore, the relation \AGn /Gn\ = 5f gives the upper limit for 5f- In 
turn, the equations (|5.ip and (|5.2p show that dp ~ l'^- Therefore, from these equations 
we can get the upper limit for /, substituting there for definiteness values for the radii 
of the spheres and the separation between them from the Moscow experiment. Thus, 
for the Washington and Ziirich experiments, in the case of the approximate formula 
(|5.ip we get respectively 9.067 /im and 10.527 /im and in the case of the exact formula 
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(|5.2p we obtain respectively 9.070 /xm and 10.531 /xm. 

Of course, we get rather rough estimates for the upper limit of I. Anyway, we 
think that it gives more or less correct value of the order of magnitude of I in the 
Randall-Sundrum model with one brane: I < 10 fim. Figure 2 shows that for such 
values of I the difference between the approximate and exact formulas is negligible. It 
is worth noting that close constraints were found in the table-top inverse square law 
experiments ^13) and from astrophysical observations |14]-|18). 

6. Conclusion 

In our paper we have considered the one-brane Randall-Sundrum model. In the 
weak-field limit, we obtained the approximate and exact expressions for gravitational 
potentials and accelerations of test bodies in these potentials for different geometrical 
configurations. Some of these approximate formulas were already known (see, e.g., 
[3 13)' exact ones were found for the first time. We applied these equations 

for calculation of the gravitational interaction between two spherical shells of finite 
thickness that can be easily reduced to the case of spheres. Then, we found the 
approximate and exact expressions for the relative force corrections to the Newton's 
gravitational force between two massive spheres. It is clear that the use of the 
approximate solution makes the calculations much easier. But we should analyze 
the difference between the approximate and exact solutions to find out where the 
application of the approximate solution is appropriate. We found that the difference 
between relative force corrections for the approximate and exact cases increases with 
the parameter I (for the fixed distance r between centers of the spheres) . On the other 
hand, this difference increases with decreasing of the distance between the centers of 
the spheres (for the fixed curvature scale parameter I). Using the results of the table- 
top Cavendish-type experiments measuring the Newton's gravitational constant, from 
the equations for the relative force corrections we got the upper limit for the curvature 
scale parameter I < 10 /im in the Randall-Sundrum model. For these values of I, the 
difference between the approximate and exact solutions is negligible. 
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